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Interpretation of radio continuum and molecular line 
observations of Sgr B2: free- free and synchrotron emission, 
and implications for cosmic rays 
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ABSTRACT 

Recent ammonia (1,1) inversion line data on the Galactic star forming region Sgr B2 
show that the column density is consistent with a radial Gaussian density profile with 
a standard deviation of 2.75 pc. Deriving a formula for the virial mass of spherical 
Gaussian clouds, we obtain Mvir = 1-9 x 10^ M0 for Sgr B2. For this matter dis- 
tribution, a reasonable magnetic field and an impinging flux of cosmic rays of solar 
neighbourhood intensity, we predict the expected synchrotron emission from the Sgr 
B2 giant molecular cloud due to secondary electrons and positrons resulting from cos- 
mic ray interactions, including effects of losses due to pion production collisions during 
diffusive propagation into the cloud complex. 

We assemble radio continuum data at frequencies between 330 MHz and 230 
GHz. From the spectral energy distribution the emission appears to be thermal at 
all frequencies. Before using these data to constrain the predicted synchrotron flux, 
we first model the spectrum as free-free emission from the known ultra compact Hii 
regions plus emission from an envelope or wind with a radial density gradient, and 
obtain an excellent fit. We thus find the spectrum at all frequencies to be dominated 
by thermal emission, and this severely constrains the possible synchrotron emission by 
secondary electrons to quite low flux levels. The absence of a significant contribution 
by secondary electrons is almost certainly due to multi-GeV energy cosmic rays being 
unable to penetrate far into giant molecular clouds. This would also explain why 
100 MeV-GeV gamma-rays (from neutral pion decay or bremsstrahlung by secondary 
electrons) were not observed from Sgr B2 by EGRET, while TeV energy gamma- 
rays were observed, being produced by higher energy cosmic rays which more readily 
penetrate giant molecular clouds. 

Key words: cosmic rays, HII regions, radiation mechanisms: non-thermal, ISM: 
individual: Sgr B2 Giant Molecular Cloud, radio lines: ISM, radio continuum: ISM 



1 INTRODUCTION 

Molecular clouds have long been studied as laboratories for 
star formation. This has led to a wealth of information about 
the physical characteristics of clouds and their chemical 
makeup, usually obtained by observing emission/absorption 
lines of molecules such as CO, OH, etc., which reveal the 
presence of molecular gas. Most molecular cloud emission is 
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thermal, from the Hll regions and/or dust emission. Cosmic 
rays play an important role in molecular cloud evolution by 
partially ionizing even the cold molecular gas, thereby af- 
fecting, through ambipolar diffusion, the dynamics of cloud 
collapse by coupling the magnetic field to the partially ion- 
ized gas. This in turn could result in amplification of the 
ambient magnetic field du ring cloud coll apse, and give rise 
to the correlation found bv lCrutcheiJ (|l999i ) between average 
density of molecular gas of molecular clouds and their line- 
of-sight magnetic field s. From V LA observations of Zeeman 
splitting of the HI line ICrutcheiJ (|l99& l found a line of sight 
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magnetic field for Sgr B2 of Blos ~ 0.5 mG to apply to the 
outer envelope of the cloud complex, and this would suggest 
that the magnetic field amplitude could be even hi gher in the 
inner parts of the complex. In fact, both iLis fc Goldsmitbl 
l|l98g ) and ICrutcherl (| 19961 ') actually countenance average 
magnetic field strengths as high as ~2 mG for the Sgr B2 
cloud on the basis of the theoretical prejudice that the cloud 
be magnetically supported against gravitational collapse; we 
cannot exclude that such field strengths may actually apply 
on large scales in the complex. 

Recent stud ies of the ionization rate by 

Ivan der Tak et al.l (|2006l ) show that the cosmic ray 
ionization rate of dense molecular clouds in the Galactic 
Centre (GC) region may be as much as a factor ~10 higher 
than in molecular clouds in the solar neighbourhood. Other 
evidence for a higher cosmic ray density in the GC region 
may come from the observation of 6.4 keV iron line emis- 
sion. Assuming that low energy cosmic rays are resp o nsible 
for heating the molecular gas, lYusef-Zadeh et al.l (|2007l ) 
estimated the energy densities of cosmic rays in GC molec- 
ular clouds to range from 19 eV cm~^ to 6 x 10* eV cm"'', 
with 51 eV cm~^ for Sgr B2. This is much higher than the 
energy density of cosmic rays in the solar neighbourhood 
which is ~1 eV cm"''. Since ionisation is most effective 
for low-energy cosmic ray nuclei, and because cosmic ray 
energy spectra typically have an inverse power-law form, 
it is the cosmic rays with kinetic energies much less than 
1 GeV/nucleon that are mainly responsible. Hence the 
enhanced ionization in the GC region could be due to an 
overall enhancement of the cosmic ray fiux there or to an 
addit ional, steep, low-energy component (jCrocker et al.l 
l200'<t ). 

High densities of low-energy hadronic cosmic ray nuclei 
within regions of dense gas will result in enhanced ~100 MeV 
gamma-ray emission through enhanced pion production fol- 



lowed by tt" -^ 77 decay and tv^ 



e decay followed 



by secondary electron bremsstrahlung. The High Energy 
Stereoscopic System ( HESS) recently compl eted a survey 
of the Galactic Centre IjAharonian et al.ll20 061. and found a 
broad scale correlation between the 7-ray emission and col- 
umn density of the molecular gas, but required the GC re- 
gion cosmic ray fiux at ~10 TeV to be about 10 times higher 
than that observed in the solar neighbourhood. Such an en- 
hancement in the GC cosmic ray flux at GeV energies was 
not inferred by EGRET observations of 100 MeV gamma- 
rays from the cen tral region of the Galaxy, indeed, i t was ex- 
plicitly noted bv iMaver-Hasselwander et al.l l| 19981 ') that no 
localized excess associated with the Sgr B complex was de- 
tected by EGRET excluding the possibility of a significantly 
enhanced CR density in these clouds - in the appropriate 
energy range, of course. One possible e xplanation of this, 
and t he high ionization rate inferr ed bv Ivan der Tak et al.l 
l|2006t ). that has been suggested bv lCrocker et al.l (120071 ') is 
that in the Galactic Centre region there may be a steep 
component i?~^'^ with a normalization at ~10 GeV energies 
comparable to that locally in order to explain the EGRET 
result, and an even steeper lower energy component to ex- 
plain the high ionization rate, and finally a flat E~^'^ com- 
ponent negligible at GeV energies but giving a tenfold in- 
crease at ~10 TeV to explain the HE SS data. Such a fla t 
component may arise as suggested by ICheng et al.l (|2007l ) 
if periodic acceleration at Sgr A* occurs when stars are 



tidally disrupted at a rate of 10~^ year~^, and diffuse at 
a distance of ~500 pc before pp losses steepen the spectrum 
on a timescale of ~10'' years. However, to also explain the 
broadband radio to gamma ray spectral energy distribution 
of the Sgr B region, with hadronic models they need a rather 
strong average magnetic field, viz., 2.2-3.7 mG. 

The same interactions of cosmic ray nuclei within re- 
gions of dense gas which may lead to enhanced gamma-ray 
production (at least at TeV energies) should produce co- 
pious secondary electrons and positrons which may in turn 
produce synchrotron emission in CMC magnetic fields which 
are observed to be higher than elsewhere in the interstellar 
medium. This possibility was the motivation for the present 
work as well as recent observations at 1.4 GHz and 2.4 GHz 
of the Sgr B2 CMC l| Jones et al.ll2008al ). and of the dense 
cold starle ss cores G333. 125-0. 562 and IRAS 15596-5301 
Jjones et a l. 2008b). The dense cores were chosen because, 
unlike Sgr B2, they are well away from the central region of 
the Galaxy and would have a cosmic ray environment ex- 
pected to be similar to that of the Solar region. They are 
of much lower mass than the Sgr B2 CMC, and their mag- 
netic fields are unknown. The non-detection of these dense 
cold starless cores in non-thermal emission was used to place 
upper limits on the magnetic fields of both of ~0.5 mG. 

Here, we investigate whether one should expect to see 
radio synchrotron emission by secondary e from giant 
molecular cloud (CMC) complexes. "We shall compare our 
predictions of the e xpected synchrotro n emission with our 
recent observations l|jones et al.ll2008ah to draw conclusions 
about the cosmic ray environment around and within the 
Sgr B2 CMC. 'We chose the Sgr B2 CMC for this study 
because of its large mass, its location in the central region 
of the Galaxy where the cosmic ray density may be higher 
than that in the solar neighbourhood, and its high magnetic 
field. This was in spite of being aware of its complicated na- 
ture, and the likely difficulty in disentangling non-thermal 
from thermal emission in this source - our search for any 
comparable molecular cloud with no star formation in the 
central region of the Galaxy was unsuccessful. 



2 SGR B2 CLOUD COMPLEX: MASS AND 
DENSITY 

The combination of magnetic fields and secondary e (and 
also primary e~) will lead to the emission of synchrotron ra- 
diation from molecular clouds which, because of its relatively 
steep spectrum, may show up at frequencies below which 
the free-free emission from Hll regions turns down after be- 
coming optically thick. The observed flux of cosmic ray e 
contains at least ^15% p ositrons at 10 GeV l|Grimani et al.l 
|2002| : iBeattv. et al.l [200J) . Given that secondary electrons 
and positrons would be produced in situ inside molecular 
clouds by cosmic ray nuclei, they should be sites of copious 
secondary e production. Since the production of these sec- 
ondaries is proportional to the product of the matter and 
cosmic ray densities within the clouds, there should be an 
appreciable flux of synchrotron radiation at low frequencies 
from molecular clouds due to secondary electrons, provided 
cosmic rays can penetrate the clouds. 

Sgr B2 is one of the largest and mos t complex molecu - 
lar cloud/Hll regions in the Galaxy - see iLang et al.l (|2008l ) 
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for a discussion and review of the continuum emission mea- 
surements. It lies near the Galactic Centre, and we assume 
it to be ~8.5 kpc from Earth, at a projected distance of 
100 pc from the Galactic Centre. Sg r B2 comprises at least 
four components (IGordon et al.lll993l ). These are three dense 
cores Sgr B2(N), (M) and (S) and a less dense outer envelope 
(OE). The dense cores are sites of massive star formation 
and have Hll regions, ultr a-compact Hll regions (UCHII) - 
ICaume fc Claussenl (|199G| ) have found more than 60 UCHII 
sources in Sgr B2(M) alone. X-ray sources associated with 
Hll regions. X-ray sources with no ra dio or IR counterparts 
(jTakagi. Murakami fc Kovamal I2OO2I ), dense cores, em bed- 
ded protostars and molecular masers (|Goicoecheall2004l ) are 
also found. The cores are small (~0.5 pc), warm («45- 
60K), light (10^-10^ Mq), dense (10'^"'' cm"^), and cor- 
respond to ~5% of the cloud mass. On the other hand, 
the envelope is cool (~10 K), massive (7.6 x IO^Mq), and 
less dense (lO^cm""^). It is thought that at wavelengths 
A >3mm (100 GHz), free-free emission dominates, whilst at 
shorter wavelengths thermal emission from dust dominates 
l|Gordon et al.lll993D . 

There have also been radio continuum, X-ray, and re- 
cently, as we mentioned earlier, 7-ray observations of the 
Galactic Centre. The large-scale diffuse radio emission from 
the GC regio n has been obse r ved at 330MHz extensively us- 
ing the VLA (jLa Rosa et al.l (|2005|) and references therein). 
It has also been obs erved in hard X-rays by INTEGRAL 
iNeronov et al.ll2005h . There have also been Chandra and 
ASCA X-ray observations of the Galactic Centre, where the 
authors argue for a reflection nebula of Compton scattered 
X-ray s from the Galactic Ce n tre source Sgr A* at an earlier 
time IIMurakami et al.| 200d : Takagi, Murakami fc Kovamal 



l2002l : iMurakami et al.ll200li:iFrver et al.ll2006i r 

Recently, lOtt et al.l ( 20061 ) have observed the ammonia 
(1,1) inversion line over the Sgr B2 complex and the result- 
ing zeroth moment map (image of the intensity integrated 
over the line) is shown in Fig. [T] The data, obtained with 
the Australia Telescope Compact Array (ATCA), show the 
Sgr B2 parent molecular cloud in sharp contrast against the 
surrounding molecular material. This is mainly due to the 
property of an interferometer to filter out very extended 
structures. Excluding the absorption of the ammonia (1,1) 
line by the prominent Hll regions Sgr B2(M) and Sgr B2(N), 
we find that the intensity integrated over the line varies with 
position in such a way that the column density profile is 
consistent with a two-dimensional Gaussian with a stan- 
dard deviation of (t=2.75±0.1 pc (assuming a distance of 
8.5 kpc) centred midway between Sgr B2(N) and Sgr B2(M) 
(see Fig. [21). 

For optically thin emission, the intensity integrated over 
the ammonia (1,1) line is proportional to the column density 
of ammonia provided the temperature is constant, and so to 
that of molecular hydrogen 

/ Iijdiy fx 7Vnh3 oc 7Vh2 , 

v'(l,l)linc 

where the final proportionality also assumes the fractional 
NH3 abundance is constant. For the assumptions above, this 
implies that the volume density profile for H2 must be a ra- 
dial Gaussian density profile, too. For the case of a cloud 
with spherical symmetry, if the column density is a two- 
dimensional Gaussian surface density, then the volume den- 
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Figure 1. Image and contours showing the zeroth moment of the 
(1,1) line of NH3 for the region around Sgr B2. J2000 coordinates 
are used. Contours from 10% peak intensity until 90% peak in- 
tensity in increments of 10%. Note that around the Sgr B2(M) 
and Sgr B(N) Hll regions, the NH3 line emission is strongly at- 
tenuated due to thermal bremsstrahlung absorption. The beam is 
located in the lower left-hand corner of the image, and is 26 X 17" 
at a position angle of —70° . The intensity scale, located on the 
right of the image is from -96 to 1280 K km s~^ . 



sity must be described by a three-dimensional Gaussian with 
the same standard deviation. 



nH2{'r) = 



Mb 



^(x^+y^+z^)/{2a^) 



2mH (^/2^CT)3 

in which case the column density at impact parameter b from 
the cloud centre is 

2mH 2Tva^ 

Of course, the column density averages over density varia- 
tions along the line of sight, and so the smooth radial Gaus- 
sian density profile will be an approximation to the true 
density distribution. Indeed, the cloud structure typically is 
fractal and the mass distribution follows a power law, and 
the column density contains contributions from a good num- 
ber of individual cloudlets. 



2.1 Virial mass of a Gaussian spherical cloud 

Here we derive, for the first time, the virial mass of a cloud 
complex with a radial Gaussian density profile. If a cloud is 
thermally supported, its kinetic energy is 



K=l 



M 



-kT 



2 {/i)m„ 

where {/i) is the mean atomic mass, m„ is the unified atomic 
mass unit, k Boltzmann constant and T the temperature. 
The mass inside radius r of a Gaussian spherical cloud is 



M{<r) 



M- 



r^/2CT^ 



1/2 -X 



dx 



M^r{3/2,r^/2a^). 



where T{a,x) is the incomplete Gamma function. The grav- 
itational potential energy of a Gaussian spherical cloud is 
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Figure 2. Azimuthally averaged intensity of the (1,1) inversion 
line of NH3 vs. impact parameter from the centre of Sgr B2 as- 
suming a distance of 8.5 kpc. The thick dashed line is a Gaussian 
fit to the data between 2 pc and 7 pc with (t=2.75 pc; at impact 
parameters less than 2 pc the data are affected by absorption ef- 
fects against SgrB2 (N) and (M), and beyond 7 pc the intensity 
may be emission from a southern cloud possibly unrelated to the 
Sgr B2 cloud (cf. Fig.[T]l. This fit provides the extrapolation of 
the profile toward the centre; the thin dotted lines give fits hav- 
ing (T=2.65 and 2.85 pc and are shown to give an indication of the 
uncertainty in a. The right hand axis shows the inferred density 
as a function of radius for the cr=2.75 pc fit. 
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Figure 3. Sketch of the morphology of Sgr B2 showing the lo- 
cations of the prominent Hll regions, the 200 K km s~^ NH3 
contour, and the locations of the main Sgr B2 cloud complex 
(shaded regions have radius equal to one and two standard devi- 
ations of the the assumed radial Gaussian density profile). The 
strong southern non-thermal source indicated is excluded from 
the present analysis. 



virial mass is given by 



then 
U = 



GMAr(3/2,rV2a^)^^^^^dr 
s/TT r 





Usually the temperature is obtained from observed thermal 
Doppler broadening of a narrow line of some element or 
molecule with atomic mass /i. Then, if the emission is op- 
tically thin, the line has a Gaussian profile with standard 
deviation (measured in m/s) of 



Cv ~ \ 


1 kT 
j fJ-rriu 


giving 




kT 


/lal. 



From the virial theorem, K = —^U, we obtain 

^^'' = ^(;i)"- 

Putting this in practical units, we obtain 






1 pc 



vfwhm 
1km 



s 



This is a factor 2.1 higher than the usual formula for a uni- 
form density sphere of radius R — a. 

If the cloud is supported solely by turbulent motion, as 
is certainly the case for GMCs, the line width is determined 
by the RMS turbulent velocity rather than the thermal RMS 
speed of the molecular species being observed, and then the 



Mvir 



= 444 



a \ / vfwhm 
1 pc 



Ikms-i 



(1) 



2.2 Mass of Sgr B2 cloud complex 

The velocity FWHM of the ammonia (1,1) Hne observations 
of Sgr B2 is 39.7 km s"'^, implying the cloud is supported 
by turbulence rather than being thermally supported. The 
virial mass of Sgr B2 based on Eq. [T]and a=2.75±0.1 pc is 
Mvir = (1.9±0.1) X lO'' M0. Since the virial mass of a cloud 
complex with a radial Gaussian density profile is a factor 2.1 
higher than that of a uniform density sphere, we suggest that 
the typical uncertainty in mass determinations using the 
virial theorem arising from unc ertainty in cloud structure 
could be as large as a factor ~2. ISatd (|200Cl ) gives a mass of 
(1-2) xlO*^ M0 for Sgr B2(M) assuming a radius of 1.5 pc 
and in the HESS paper on the GC region lAharonian et al.l 
(2006) give a total mass of (6-15) x 10*^ Mg for a 0.5° x 0.5° 
= 75 PCX 75 p c region surrounding Sgr B2 based on CS data 
(TsuboilllQQa ). 

Taking the virial mass of (1.9 ± 0.1) x lO'^ Mq to be 
the mass of molecular gas (assumed to be almost entirely 
H2) in the Sgr B2 cloud complex, and a radial Gaussian 
density profile with a standard deviation of cr=2.75±0.1 pc, 
the maximum column density is ^h2 = (2.5 ±0.1) x lO'^* 
cm"^. The density is nj/^ = (1.2 ± 0.1) x 10^ cm~^ at the 
centre of the Sgr B2 complex, and decreases to 10 cm~^ at 
a radius of ~12 pc which we shall consider to be its outer 
radius. The H2 number density (dotted curve) may be read 
off the right-hand axis in Figure [5] A sketch of the Sgr B2 
region showing the ammonia (1,1) hne 10% contour level of 
the zeroth moment map, the (N), (S) and (M) Hll regions. 
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and the size of inferred Gaussian cloud complex is given in 
Figure [1 



3 COSMIC RAY SECONDARY ELECTRON 
PRODUCTION 

The Galactic synchrotron emission is due to accelerated (pri- 
mary) cosmic ray electrons, and to electrons and positrons 
produced as secondaries. The production rate of secondary 
electrons and positrons depends only on the spectrum and 
intensity of cosmic ray nuclei, and the density of the in- 
terstellar matter. We use the production rate of electrons 
and positrons qm(E), per solar mass of interstellar mat- 
ter per unit energy (Mq^ GeV~^ s~^), for the cosmic ray 
spectrum and composition observed above the Eart h's at - 
mosphere based on figure 4 of Moskalonko & Strong ([1993) ■ 

The production of large numbers of positrons can be 
traced through their 511 keV annihilation line, and this has 
been observed by INTEGRAL from with in 8° of the Galac- 
tic Centre region (Wcidenspointncr 2 0061 ) . To explain these 
observations, Beacom & Yiiksel (200^) show that the pro- 
duction of positrons is ~ 10^° per year in the central region 
of the Galaxy, and that they must be injected at energies 
below ~3 MeV to avoid excessive gamma ray emission at 
higher energies (excluding the possibility that they are sup- 
plied by the pp — > tt"*" — > e"*" chain). 

The production spectrum of secondary cosmic ray e 
has a gradual cut-off below ~0.3 GeV due to threshold 
for TT in p — p collisions. Nevertheless we should check 
that their production rate in molecular clouds does not 



exceed the stringent constraint on production of 



10=' 



positrons per year. Assuming the cosmic ray spectrum in 
the Galactic Centre region has the same shape there as 
locally, but is enhanced by a factor fcR, we find the to- 
tal production rate of secondary e^ in Sgr B2 alone to be 
1.9 X 10** MQfcRfqm{E)dE = 7.8 x 10*Vch e+/year. 
This certainly does not exceed the limit for fcR in the 
range 1-10 which seems reasonable based on EGRET and 
HESS gamma-ray data, and taking account of the likeli- 
hood that cosmic rays may not fully penetrate clouds the 
production rate is probably less. Indeed, one would require 
~ 3 X 10^^ Mq/ fcR of interstellar gas within 1 kpc of the 
Galactic Centre to exceed the observed e^ production rate. 
We shall delay until later in this section a discussion 
of the complex problem of penetration of Galactic cosmic 
ray nuclei and electrons into molecular clouds. For the mo- 
ment, then, we shall assume that Galactic cosmic rays freely 
penetrate the cloud and that their spectrum inside the cloud 
complex is the same as that in the solar neighbourhood. The 
production spectrum of electrons and positrons per unit vol- 
ume per unit energy (e cm~^ GeV~^s~^) at position r in 
the Galactic Centre region is then the product of the density 
of interstellar gas at position r multiplied by the production 
rate per unit mass 



q{E,r) = fcRnH2(r)qrniE)2mH/MQ. 



(2) 



For moderate molecular cloud densities uh > W^ cm^'^ and 
magnetic fields B > 10~^ G, the relatively short energy 
loss times appear to justify neglecting diffusive transport of 
electrons - we shall discuss this point in detail in a later 
section. Then one readily obtains, by numerical integration. 



the ambient number density of electrons and positrons, per 
unit energy, n {E, r)(e cm~"^ GeV~^), at various positions 
r within the molecular cloud complex: 

rq{E,r)dE 

where dE/dt is the total rate of energy loss of electrons 
at energy E due to ionization, bremsstrahlung and syn- 
chrotron emission (because of the energies involved, we 
neglect positron annihilation, and assume electrons and 
positrons suffer identical energy losses). Electrons lose en- 
ergy by ionization losses in neutral molecular hydrogen at a 
rate (in GeVs~^) of 

^ =5.5xl0-^^f^^ 

dt ioniz V 1 Cm~ 

and by bremsstrahlung at a rate (in GeVs"^) of 

^b.on.ss " ^-^ "" ^°''' (iGiv) "" (l^^j ■ 

The synchrotron energy loss rate (in GeVs"^) is 



3^ X (ln7 + 6.85) 



^ = 1.0 X 10-^^ X , , 

dt synch V 1 gauss 



X 7 



where B± is the component of magnetic field perpendic- 
ular to the electron's direction. For an isotropic electron 
population the solid-angle average is {B±} — 7r_B/4, and, 
assuming the magnetic field can be in any direction with 
respect to the line of sight, an appropriate value for the 
line-of-sight component of magnetic field obtained from Zee- 
man splitting Blos would be {B^os) = B/2. Hence taking 
B± = iiBlosI'2 is reasonable. Given the observed value is 
Bios=0.5 mG (|CrutcheiJll99^ ). we adopt Bx=0.8 mC. 

The synchrotron emission is calculated using 
standard formulae in synchrotron radiation theory 
iRvbicki fc Lightman|[l979h 



J^ = 



B, 



VSe^ 



i-KrUeC^ \ 1 gauss 



erg cm '^ s ^ sr ^ Hz 



Uc = 4.19 X 10^{E/meC^f 



Bi 



1 gauss 



F{v/uc)n{E,r)dE 



Hz, 



4.8 X 10" 



m^c- = 8.18 X 10~^ 



erg. 



F{x) 



K.{i)di. 



and K^ (x) is the modified Bessel function of order 5/3. 

The Razin efi'ect reduces non-thermal emission at low 
frequencies by suppression of synchrotron emission a.t u < 
7e!^p where Vp is the plasma frequency and 7e is the Lorentz 
factor of the radiating electrons, and the effect is negligi- 
ble where v 3> 20ne/B where Tie is the number density of 
free electrons (cm~ ) and B is in Gauss. For the Sgr B2 
cloud assuming i3x=0.8 mG and our lowest frequency of 
interest being 330 MHz, the Razin effect will be small if 
Tie ^ 10^ cm""^. Given that in a molecular cloud Ue ^ nH2 
and for the Sgr B2 main complex we estimate the central 
density to be nH2 = 1-2 x 10^ cm~^, we can safely neglect 
the Razin effect in the present work. 

For a perpendicular component of magnetic field of 
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Figure 4. Synchrotron emission coefficient of secondary etectrons 
produced by cosmic ray interactions for densities nu^ = lO" (bot- 
tom curve), 10^, 10^, . . . 10® cm""^ (top curve). 



0.8 mG, as may be appropriate for Sgr B2, and the so- 
lar neighbourhood cosmic ray spectrum, we find the spe- 
cific emission coefficient due to synchrotron emission by sec- 
ondary electrons as shown in Fig.Ufor various gas densities. 
Note that for higher density regions, electron energy loss by 
bremsstrahlung dominates over synchrotron losses reducing 
the synchrotron power per unit mass relative to lower den- 
sity regions. As can be seen, this has the effect that at low 
frequencies the synchrotron emission coefficient becomes al- 
most independent of density for nu^ > 10^ cm~'^ and mG 
fields. 



3.1 Penetration of Cosmic Ray Nuclei into the 
Sgr B2 Cloud Complex 

Work on this subject has been motivated mainly by gamma- 
ray observations, particularly of the central region of the 
Galaxy. An important contribution to the Galactic gamma- 
ray intensity comes from interactions of cosmic ray nuclei 
through pion production and subsequent decay tt'' —* 77, 
and TV —* fi -^ e followed by bremsstrahlung or inverse 
Compton. Of course primary accelerated electrons are also 
important for the latter two processes. Put simply, if cosmic 
rays can freely enter molecular clouds then the gamma-ray 
flux will be higher than if they cannot. In the present work 
we are interested in synchrotron radiation by the same sec- 
ondary e . Another motivation has been to estimate the 
ionization rate due to cosmic rays. Again, this depends cru- 
cially on the extent of penetration of cosmic rays responsible 
for ionization, mainly those of lower energy. 

The nature and extent of penetration of cosmic rays 
into molecular clo uds has is not yet fully understood, and 
has a long history. ISkilling fc Strong (|l976l ) concluded that 
the very low energy cosmic rays mainly responsible for 
ionization of clo ud material are e f ficient ly excluded from 
clouds, whereas ICesarskv fc Volkl (11978! ) concluded that 
molecular clouds are p ervaded by inter-cloud cosmic rays. 
iDogel' fc Sharovl (|l990l ) considered acceleration of charged 
particles by turbulence in giant molecular clouds, and sug- 
gested that this mechanism may explain the unidentified 
gamma-ray sources discovered by COS-B, and estimated the 



synchrotron radio emission of accelerated primary and sec- 
ondary electrons in molecular clouds in this scenario. 

Certainly, at multi-TeV energies there should be no 
problem in cloud penetration, and this was recently con- 
firmed by the excellent correlati on between TeV gamm a-ray 
intensity as measured by HESS (jAharonian et al.ll2006l ) and 
the column density of molecular gas for the Galactic Centre 
region. At lower energies, the interpretation of the 100 MeV- 
GeV energy gamma ray intens ity measured by EGRET 
(jMaver-Hasselwander et al.lll998l ') toward the Galactic Cen- 
tre region is ambiguous and we await with great interest the 
higher resolutio n data from GLAST due for launch in 2008. 
In this context, ICabici et al.l (|2007l ) have recently investi- 
gated the penetration of cosmic rays into molecular clouds to 
understand the importance of this for gamma ray emission 
at GeV and TeV energies. They considered proton-proton 
collision losses as cosmic rays diffuse into a cloud. Taking 
a typical cloud to have n//2~300 cm~^, _B=10 /iG and ra- 
dius 20 pc, they found that, for a diffusion coefficient based 
on that which seems to apply to cosmic rays throughout 
the Galaxy (as determined by secondary to primary compo- 
sition measurements), cosmic rays would freely penetrate. 
They also found, however, that if the diffusion coefficient 
inside the cloud is smaller, say 0.01 of the average Galactic 
one, that exclusion becomes relevant for 10-100 GeV cosmic 
ray nuclei resulting in suppression of GeV gamma-ray emis- 
sion. Given that the Sgr B2 complex has a muc h higher den- 
sity an d magnetic field than that modelled bv lCabici et al.l 
(|2007|) , it is clearly necessary to determine the extent of sup- 
pression of e production. 

Cosmic ray protons and nuclei produced outside the 
cloud will penetrate the cloud by diffusion and advection, 
and lose typically half their energy in pp collisions on a time 
scale tpp ~ 5 X lO'^n^^ yr, where 71^2 is in cm~^. The ad- 
vection timescale is fadv = Rcioud/crv which for i?cioud=12 
and ay— 39.7 km s~^ is tadv — 3 x 10® years. 

The diffusion timescale is 



pc 



Uifi{E) 



R; 



2 

cloud 



(3) 



2D{E) 

where D{E) is the diffusion coefficient, which depends on 
the ambient magnetic field and the spectrum of turbulence. 
The minimum diffusion coefficient for a completely tangled 
magnetic field is the so-called "Bohm diffusion coefficient" 
which, for relativistic protons, is -Dmiii(-E) — ^rg{E)c ex E 
where Vg ~ 10~^EgcvB^q pc is the gyroradius, BmG is the 
magnetic field in milligauss and -EceV is the proton energy in 
GeV. Models of cosmic ray propagation in the Galaxy which 
are consistent with the observed relative abundance of "pri- 
mary" cosmic ray nuclei (e.g. Carbon) and "secondary" cos- 
mic ray nuclei (e.g. Boron) - the latter produced by spalla- 
tion of primary cosmic ray nuclei - suggest that D{E) oc E" 
where a ~ 0.3 — 0.7. If a Kolmogorov spectrum of turbu- 
lence is present, then one would expect a— 1/3, but in the 
presence of strong magnetic fields a Kraichnan spect rum 
may give rise to a= 1/2. Following iGabici et al.l ((20071), we 
adopt a diffusion coefficient 



D{E) = 3 X 10^ X 



E/{1 GeV) 
S/(3 /iG) 



(4) 



where x =^ 1 is a factor to account for the possible suppression 
(slowing) of diffusive transport. 
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Adopting a magnetic field of ~0.8 mG, tfie typical en- 
ergies of electrons or positrons responsible for synchrotron 
emission at 330 MHz - 1 GHz are ~0.3-0.4 GeV. Taking 
the primary proton energy to be ~10 times higher, and set- 
ting the outer boundary of the Sgr B2 cloud to be ~12 pc, 
this gives a diffusion timescale of ~ 10*x~^ yr- Comparing 
this diffusion time-scale with the loss timescale of only 500 
years, for pp collisions in a central density of nH2 ~ 10^ 
cm~^, shows that, for the 0.8 mG magnetic field, penetra- 
tion to the centre of the GMC complex is very improbable 
below ~3 GeV energies. As we have already noted, because 
bremsstrahlung losses dominate in dense regions much of the 
synchrotron emission is expected to come from outer regions 
of Sgr B2, and it is penetration to these outer regions that 
matters most. For these regions the distance is obviously 
smaller and the density lower, suggesting that penetration 
of the outer region of the GMC complex is less of a problem. 

We can make a more quantitative approximation of 
cosmic ray nucleus penetration by analogy with scatter- 
ing ('s') and absorption ('a') of radiation and we define an 
effective optical thickness t-^ — Taija + Ts) analogous to 
that used when cons i dering radiative diffusion - see, e.g., 
iRvbicki fc LightmanI (119791 ). In our case, for penetration 
from an outer boundary 7? to distance r from the centre 
we have 



Ta{r) 



0.5[2nH2{r')]o-ppdr' 



(5) 



where app{E) «35 mb above threshold, and the 0.5 factor is 
approximately the mean inelasticity (fractional energy lost) 
in pp collisions, and 



s{E, 



-dr 



(6) 



3D{E, 

since for isotropic diffusion, the effective mean free path is 
3D/c. Then the cosmic ray intensity at radius r is Icr ~ 

The mean primary proton energy for a given secondary 
electron energy is given by 

jEpnp{Ep)Y{Ee;Ep)dEp 



(Ep) 



Jnp(Ep)Y{E,;Ep)dEp 



where np{Ep)dEp (cm~'^ GeV~^) is the number density and 
of cosmic ray protons with energy Ep to {Ep + dEp), and 
Y{Ee; Ep)dE^± (g~^ cm^) are particle yields giving the rate 
of production of secondary electrons and positrons with en- 
ergy Ee to (-Ee + d-Be) per unit pathlength (in g cm~'^) by 
a single cosmic proton of energy Ep. These particle yields 
can be obtained from accelerator data on charged pion pro- 
duction in pp collisions taking account of tt -^ /i — > e decays 
(we use data kindly provided by T. Stanev, private commu- 
nication). We find that for the local cosmic ray spectrum 
the mean primary proton total energy is 



{Ep) « 0.0157e + 227; 



GeV. 



(7) 



So for observations made at frequency i/, the appropriate 
cosmic ray proton energy to use is determined by assuming 
electrons radiate at the critical frequency, i.e. putting 

1/2 



7e = 



,4.19 X 106 Bi 
into Eq. [T] where v is in Hz, B± 



(8) 



l:^ 






1.0 








....-■■.■;.--^-? :f 


0.8 


thin curves: 5^=0.8 mG 
thick curves: 5i=0.1 mG/ 


■ ■ ' /' / / . 


0.6 


■ x»i 


■■' / '' ' / 


0.4 


X=l 

x=0-oi / . / 




0.2 


/ 


/ /' 


no 







a 4 6 8 10 12 

r (pc) 

Figure 5. The cosmic ray penetration factor e"'^*!^?''''''''''! ap- 
propriate to ;/=0.3 GHz is plotted against radius for the diffusion 
coefficient defined by Eq|4]with magnetic field and x '^^ labeled 
(x 2> 1 corresponds to unimpeded penetration). 

The penetration factor e~^*'^''' appropriate to an ob- 
serving frequency of 0.3 GHz is plotted against radius from 
the centre of the Sgr B2 GMC for various values of the diffu- 
sive transport suppression factor x in Fig.[5]for B_l=0.8 mG. 
To show how our results depend on assumed magnetic field, 
here and elsewhere in this paper we shall also show results for 
a lower magnetic field, arbitrarily chosen to be B_l=0.1 mG. 
The penetration factor for this magnetic field is also shown 
in Fig. El 

3.2 Diffusion of Secondary and Primary Cosmic 
Ray Electrons in the Sgr B2 Cloud Complex 

Before predicting the synchrotron emission from the Srg B2 
GMC, we consider the diffusion of secondary e and primary 
cosmic ray e~. Due to the high magnetic field, and high 
densities within the inner part of the GMC, electrons will 
suffer rapid energy losses there and this will limit how far 
they can propagate by diffusion. An estimate of how far an 
electron with Lorentz factor 7e can propagate by diffusion 
before losing a significant fraction of its energy is given by 
what we shall refer to as the "diffusion-loss distance" 

1/2 



^ditt{le 



D{~femeC )7e 
(d7e/dt)total 



(9) 



m gauss. 



We plot a;5iff (7e) for B±=0.8 mG and x=l in Fig-El and also 
for i3j_=0.1 mG to show the effect of a significantly lower 
magnetic field than appears to be present over the Sgr B2 
GMC. For lower x values, the diffusion-loss distance is lower 
and is obtained by multiplying by ^/x. 

For synchrotron radiation at the adopted magnetic field, 
i3_L=0.8 mG, the Lorentz factor of electrons mainly respon- 
sible for emission at 1 GHz is 7e«550. Near the centre of 
the GMC where nH2^W^ cm~ the diffusion loss distance 
is x^iff(550) ~0.2 pc, at r=6 pc where nH2~10* cm""^ we 
find a::diff(550) ~0.5 pc, at r=8 pc where ni:/2~2xl0^ cm~^ 
we find ~1.3 pc, and at r=12 pc where nij2~10 cm~^ we 
find a:|i(j(550) ~3.7 pc. From this we can draw the following 
conclusions: (a) at all distances the diffusion loss distance of 
electrons producing synchrotron radiation at 1 GHz is small 
compared to the radial coordinate and so we may safely 
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Figure 6. The diflfusion-loss distance as defined by Eq. |9]plotted 
vs. 7e for the diffusion coefliicient given by Eq 3] for two mag- 
netic fields (as labelled) and n//2=10^ cm~^ (bottom curves), 
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Figure 7. Emission coefficient at frequency i/=0.3 GHz vs. dis- 
tance from the centre of the Sgr B2 GMC for B^ as indicated 
and for diffusion as defined by Eq|4]with x as labeled. 



make the approximation that secondary electrons within the 
GMC radiate where they are produced; (b) primary cos- 
mic ray electrons from outside the GMC will not be able to 
propagate significantly towards the centre of the GMC and 
so are effectively excluded from the GMC complex. These 
conclusions are made even stronger if x<l- However, pri- 
mary electrons accelerated inside the GMC, e.g. by diffusive 
shock acceleration at supernova shocks or wind shocks, will 
produce synchrotron emission inside the GMC. 

For the case of a weaker magnetic field, e.g. 
5^=0.1 mG, which is lower than appears to be present 
in the Sgr B2 GMC but which may occur in some other 
clouds, the Lorentz factor of electrons mainly responsible 
for emission at 1 GHz is 7e~1500, and more significant pen- 
etration of primary cosmic ray electrons would take place 
unless X ^ 1- 



3.3 Predicted synchrotron emission from Sgr B2 

At frequency v, for each point within the Sgr B2 giant molec- 
ular cloud complex we determine the H2 number density to 
find the synchrotron emission coefficient corresponding to 
complete penetration of cosmic rays within the cloud com- 
plex. Multiplying this by the cosmic ray penetration factor 
and fcR we obtain the synchrotron emission coefficient j^ (r) 
taking account of cloud penetration and the possibility of 
cosmic ray enhancement in the Galactic Centre region com- 
pared to that locally. The synchrotron emission coefficient 
at 0.3 GHz is plotted as a function of distance from the cen- 
tre of the Sgr B2 GMC for Bx=0.8 mG and Bx=0.1 mC 
in Fig. [7] In both cases, for x ^ 1 the emission coefficient 
at the cloud complex centre is negligible compared to that 
near its edge. 

We obtain the intensity Iv(&), shown in Fig.[8]as a func- 
tion of angular distance 6 for two frequencies, by integrating 
through the cloud complex assuming the synchrotron emis- 
sion is optically thin. 



Depending on the diffusive transport suppression factor x. 



we may expect significant "limb brightening" of the syn- 
chrotron emission. We find that in the case of the Sgr B2 
complex most of the flux comes from within ~11 pc of its 
centre. 



4 DISCUSSION 

Radio continuum observations of the Sgr B2 GMC, including 
new measurements at 1.4 GH z and 2.4 GHz, are d iscussed in 
detail in a separate paper by I Jones et al.l l|2008ah . Since we 
want to compare the continuum emission in the Sgr B2 re- 
gion with our predictions for the region of the dense central 
region of the giant molecular cloud we use continuum fiux es- 
timates for the region 11 pc in size centred midway between 
Sgr B2 (M) and Sgr B2 (N). Note that other papers may use 
quite different sizes for the Sgr B2 region and hence quote 
very different fluxes. In this complex region of the Galaxy, 
when assembling radio spectra from fluxes established by 
different groups, it is essential to use fluxes obtained over 
the same solid angle. 

No evidence was found for diffuse, non-thermal emis- 
sion out to ~11 pc from the centre of the Sgr B2 GMC with 
limb-brightening as predicted for synchrotron emission by 
secondary el ectrons in the pre vious section, consistent with 
the results of lLang et al.l l|2008r ) who flnd a thermal spectrum 
for Sgr B2 including its envelope. Nor was there any sug- 
gestion of polarised emission characteristic of synchrotron 
emission. However, the major Hll regions Sgr B2(M) and 
Sgr B2(N) showed up at all frequencies, and there is evi- 
dence of a strong non-thermal source to the south, which we 
have marked as "Non- Thermal Source" in Fig. [31 which is 
the subject of a separate paper ?. Figure [9] shows the spec- 
tral energy distribution for the central region of the Sgr B2 
GMC (excluding the Non-Thermal Source), and therefore 
includes the combined fluxes from the major Hll regions. 

Sgr B2 is by far the most massive molecular cloud in the 
Galaxy and it is near the Galactic Centre which is almost 
certainly a region of enhanced cosmic rays. When selecting 
Sgr B2 for our study, we knew that separating out the ther- 
mal emission would be extremely difficult. Our search for 
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Figure 8. Predicted specific intensity Iv{6) at (a) 330 MHz, and 
(b) 1.4 GHz, as a function of angular distance from the centre of 
the Sgr B2 complex for diffusion as defined by Eq |4] with x as 
labeled, a cosmic ray flux equal to that at Earth and for magnetic 
field as labeled. 



any comparable mass molecular cloud with no star forma- 
tion was unsuccessful. 

Before estimating upper limits to any synchrotron emis- 
sion from secondary electrons in the Sgr B2 GMC, we shall 
attempt to fit the observed spectral energy distribution 
solely by thermal emission processes. A large number of 
UCHII regions have been observed at high frequencies, and 
so our first step will be to extrapolate their spectra to low 
frequencies. This will be done under the assumption that 
each known UCHII region is a homogeneous sphere of ion- 
ized interstellar gas, and we shall sum the contributions from 
all known UCHII regions. Clearly this will be an approxi- 
mation as there will also be contributions from as-yet undis- 
covered UCHII regions, and because many of the UCHII 
regions will not be homogeneous, having density gradients, 
e.g. winds. We shall find that the UCHII regions account 
for ^50% of the high frequency fiux, and give a negligible 
fraction of the observed flux at low frequencies under these 
assumptions. The second step will be to fit the residual flux 
as thermal emission. The shape of the spectral energy dis- 
tribution between 330 MHz and 1.4 GHz is suggestive of 
thermal emission from a region, or regions, with a density 
gradient, and we shall model it as a free- free emission from 
one, or many identical, single temperature winds. While this 
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Figure 9. Observed fluxes summarized by Ijones et al.l ll2008al ) 
from the central region of Sgr B2 complex including the major 
Hll regions but excluding the Southern Non-Thermal Source. The 
flux from the known UCHII regions is indicated (chain curve), 
and the best fitting model of free-free emission from a constant 
temperature spherical envelope or wind with n oc r~^ is shown 
by the dashed curve, and the solid curve gives the sum of the two 
thermal components. 



is clearly unrealistic, the data available to us do not justify 
a more sophisticated treatment. 



4.1 Free- free emission from UCHII regions 

The emission at the higher frequencies (22 GHz and 43 GHz) 
is clearly thermal. Regions with very low emission measures 
and high 22 GHz and 43 GHz fluxes could potentially af- 
fect the emission at low frequencies. In order to investi- 
gate how much of the fiux at these lower frequencies could 
be attributed to UCHII regions, we modeled the emission 
from the ^ 60 known indiv i dual c omp act and UCHII region s 
reported in ICaume et all (jiggsl ) and Ide Free et all l|l998h . 
This was achieved by "bootstrapping" the fiux at the re- 
spective frequencies such that 



S^ 



k 



V 

IF) 



1 — e " 



where Vi — 22 GHz or 43 GHz, the sum over the 
~60 UCHII regions with the label (fe) relates to the fcth 
UCHII region, and the frequency dependence of the ther- 
mal bremsstrahlung absorp tion coefficient is taken from 
iRvbicki fc LightmanI l|l979 ). The summed fiux from the 
known UCHII regions has been added to Fig. [O] and it can 
be seen that the total flux from these regions accounts for 
about 50% of the 23 and 230 GHz fiux but their contribution 
below 3 GHz is negligible. 



4.2 Free-free emission from envelopes or winds 
with density gradients 

Between 330 MHz and 1.4 GHz the spectrum may be fitted 
with a single power-law Si, ~ v'^'^ characteristic of optically 
thick emission from a spherical envelope or wind with a den- 
sity gradient of the form 
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ne{r) = ni{r) = no ( — ] 



as described by IPanagia fc Fellil l|l97a) who give the ex- 
pected flux at low frequencies 



S-- = 0.611 f^^V'-I^ f-^ 

Vlcm-^y Vlpc/ VIOG 
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GHz 
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For the case of optically thin emission we can use 



4:Tvr ne{r)ni{r) dr = 4nno tq 



together with the free-free emission coefficien t and Gaunt 
factor g{T, v) from lRvbicki fc LightmanI (|l979l ) to obtain the 
flux at high frequencies where it is expected to be optically 
thin 



Sr = 1.6xl0-^f^^) -^ 
Vlcm"''/ ylpc 

/ T \-°-'' , J d 



Taking the optical depth to be r^ = ST'^jST the flux 
from the envelope or wind is Sv — 5'^'^"^'*(1 — e^'^"). 

We fit this density-gradient model to the observed 
330 MHz to 230 GHz fluxes (after subtracting the contri- 
butions of known UCHH regions as shown by chain curve in 
Fig. ^ and this is shown by the dashed curve in Fig. [^ the 
solid curve shows the sum of the two thermal components. 
For a temperature T = lO** K the best fitting parameters 
are no = 3.47 x lO'' cm~^ and tq = 4.12 x 10"^ pc. The 
high density and small size would indicate that the emis- 
sion is likely to have come from winds off, or excited by, 
young stars within the Hll regions. If the fiux is due to A'^ 
separate identical objects, their wind parameters would be 
no = 3.47 X lO'^ X iV cm"^ and ro = 4.12 x 10"^/iV pc in 
order to give the same total spectrum. 



4.3 Synchrotron emission by secondary electrons 

In Fig. [To] we re-plot the spectral energy distribution and 
show the fiux predictions for the synchrotron emission from 
secondary electrons. Given that we have found no evidence 
of synchrotron emission from secondary electrons, and that 
the observed radio continuum emission is consistent with a 
thermal origin, we shall require that any fiux of synchrotron 
emission from secondary electrons be well below the ob- 
served 330 MHz fiux. For this, we shall somewhat arbitrarily 
adopt an upper limit of 5,^ = 1 Jy at 330 MHz for any non- 
thermal component. 

Comparing the predicted synchrotron fiux in Fig. llOf al 
with our adopted limit, we find that for _B^=0.8 mG and 
fcR=l, the diflusion factor must be x <0.02. With a sig- 
nificantly lower magnetic field such as _B^=0.1 mG (see 
Fig. llOb ). even x ^ 1 (unimpeded cosmic ray penetration) 
is allowed and for this case cosmic ray enhancement in the 
Galactic Centre region at multi-GeV energies up to a fac- 
tor fcR <2.5 higher than in the Solar neighbourhood is 
not excluded. For this magnetic field, and more reasonable 
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Figure 10. Observed fluxes summarized by I Jones et al.l ll2008al ) 
from the central region of Sgr B2 complex including the major Hll 
regions but excluding the Southern Non-Thermal Source plotted 
together with fluxes predicted for (a) B± = O.SmG, (b) B± = 
O.lmG, and a cosmic ray spectrum as in the solar neighbourhood. 
Predicted fluxes are shown for diffusion as deflned by Eq|4]with 
X as indicated. The total estimated thermal flux from the major 
Hll regions is shown by the upper solid curve. 



diffusion factors, the maximum allowed cosmic ray enhance- 
ment is 3 (x=l), 3.9 (x=0.1), 7.7 (x=0.01). We emphasize 
that these are upper limits for cosmic ray enhancement and 
that there is no evidence for any cosmic ray enhancement at 
multi-GeV energies. In fact for the higher magnetic field, i.e. 
B_L=0.8 mG, which is based on Zeeman splitting observa- 
tions, the data suggest that there is no enhancement, or that 
cosmic rays at these low energies are unable to significantly 
penetrate into the Sgr B2 CMC. 

In conclusion, we have no evidence that synchrotron 
emission by electrons and positrons produced by cosmic ray 
interactions has been observed from the Sgr B2 molecular 
cloud complex. The most likely explanation for this is that, 
for reasonable diffusion models, cosmic rays with multi-GeV 
energies (that produce secondary electrons with the right 
energy to radiate at GHz frequencies in ~0.8 mG fields) 
cannot penetrate into the dense central regions of Sgr B2 
GMC where much of the potential mass of target nuclei 
is located. This exclusion is also the likely explanation for 
non-observation of the Sgr B2 GMC by EGRET because it 
is again the same multi-GeV energy protons producing pious 
in pp collisions followed by vr" -^ 77 that make an important 
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contribution to 100 MeV to multi-GeV Ramnia-rays. The 
obser vation of the Sgr B2 GMC by HESS (|Aharonian et alj 
120061 1 at TeV energies is consistent with more complete pen- 
etration of cosmic rays at higher energies into the dense cen- 
tral regions. 

In choosing giant molecular clouds in the central region 
of the Galaxy for future investigation of their synchrotron 
emission by secondary electrons, one would look for a GMC 
with with a mass of a few 10^ Mq , a lower central density 
than Sgr B2, e.g. nH2~10^ cm"'' so that low energy cosmic 
rays may more easily penetrate it, a magnetic field above 
0.1 mG and little star formation. We do not know of any, 
but such clouds may become apparent with the aid of new 
infrared surveys. 

Finally, we emphasise that as we have no independent 
knowledge of the diffusion coefficient within the Sgr B2 
GMC, i.e. X is unknown, we are unable to make a definitive 
statement about the enhancement of the low energy cosmic 
ray flux in the central region of the Galaxy relative to that 
in the Solar neighbourhood. 
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